Crossed Andreev Reflection in Diffusive Contacts 
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Crossed Andreev reflection in multiterminal structures in the diffusive regime is addressed witliin 
the quasiclassical Keldysh-Usadel formalism. The elastic cotunneling and crossed Andreev reflec- 
tion of quasiparticles give nonlocal currents and voltages (depending on the actual biasing of the 
devices) by virtue of the induced proximity effect in the normal metal electrodes. The magnitude 
of the nonlocal processes is found to scale with the square of the barrier transparency and to de- 
cay exponentially with interface spacing. Nonlocal cotunneling and crossed Andreev conductances 
are found to contribute equally to the nonlocal current, which is of relevance to the use of normal 
metal-superconducting heterostructures as sources of entanglement. 

PACS numbers: 03.67.Mn, 73.23.-b, 74.45.+C, 74.78.Na 



I. INTRODUCTION 

In standard Andreev reflection at a single normal 
metal-superconductor interface an electronlike quasipar- 
ticle in the normal metal can be transformed into a hole- 
like quasiparticle of opposite momentumii When two 
normal metal [N] electrodes or ferromagnets [F] are at- 
tached to a superconductor [S) at a distance from each 
other of the order of the coherence length, two addi- 
tional nonlocal processes are possible. During elastic co- 
tunneling (EC), an electron is transferred from one elec- 
trode to the other, while for crossed Andreev reflection 
(CAR) an electron in one of the electrodes is transformed 
into a hole in the other electrodei^i^ 

Bell-inequality experiments, quantum computation, 
and teleportation of quantum states require quantum ob- 
jects that are entangled. Cooper pairs in superconductors 
are spin singlets and are, therefore, suitable sources of 
entanglement. Crossed Andreev reflection is a promising 
possibility for the creation of locally separated entangled 
electrons iii^i^iii^ 

Nonlocal transport properties have been seen exper- 
imentally in NS^^ and FS heterostructures i^^'^^ The 
microscopic origin of the effects of EC and CAR, as well 
as possible ways in which these can be used for the pro- 
duction of locally separated entangled quasiparticles, re- 
main to be understood theoretically. 

Recently, theoretical studies of CAR have been per- 
formed for various types of j unctions . ^^'^"^i^^'^^i^^'^^i^^ 
Modeling of nonlocal effects by means of perturbation 
theory using an effective tunnel Hamiltonia n^'^i has in- 
deed provided EC and CAR, the signal being, however, 
vanishingly small (of the order of the fourth power in in- 
terface transparency). Consequent pioneering efforts to 
include disorder— and weak localization!^ have found 
enhanced effects, but still not to the level of experimental 
observations, so that the question as to the microscopic 
origin of CAR still remains open. 

In this paper, we provide insight into the underly- 
ing microscopic mechanism of CAR by studying nonlo- 
cal transport by means of quasiclassical kinetic theory. 
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FIG. 1: (a) Schematic representation of the double-barrier 
NSN structure, which is voltage biased accross the first in- 
terface. An incoming electron from the left normal metal elec- 
trode, can undergo three different processes that contribute 
to the current: Andreev refiection (AR), elastic cotunneling 
(EC) and crossed Andreev refiection (CAR). The depicted 
spectral currents that result from these processes are the main 
result of this paper, (b) Results have been obtained from a 
one-dimensional Keldysh-Usadel quasiclassical Green's func- 
tion calculation of the depicted structure. 



We present a mechanism in which CAR exists by virtue 
of the proximity induced superconducting correlations in 
the electrodes, for arbitrary barrier transparency. The 
proximity effect is the essential new ingredient in our 
model, giving a large contribution to CAR, of second or- 
der in transparency, in contrast to the nonlocal effects 
of tunnel Hamiltonian models. We show how CAR re- 
lates to the competing process of co-tunneling. For bal- 
listic transport, Andreev reflection is understood most 
straightforwardly, but we have modelled the additionally 
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challenging case of diffusive transport, as experimentally 
often is the case. Our model is of relevance to the fu- 
ture design of experiments that are based on Andreev 
entanglers. 



II. QUASICLASSICAL MODEL 

The most generic model system to study nonlocal ef- 
fects in superconducting structures is a three-terminal 
configuration consisting of a quasi-lD superconducting 
wire of length d, attached to normal reservoirs Ni^2 and a 
superconducting reservoir S which can be independently 
biased, see Fig. [T] This model is an extension to the ear- 
lier approach by Volkov et al?'^ that was used to calcu- 
late nonequilibrium transport properties of two-terminal 
N'NS contacts. The electrodes Ni^2 are weakly coupled 
to the wire, while the reservoir S is in good electrical 
contact with the wire so that their electric potentials are 
equal. Morten et alt^ addressed in their circuit model 
the role of the latter coupling strength between supercon- 
ductor S and wire. They found nonvanishing nonlocal 
effects, but only in the case of weak coupling. In the case 
of zero resistance between wire and S (strong coupling) , 
the circuit results coincide with the tunnel Hamiltonian 
results of a vanishing nonlocal signal. However, in the 
experiments often no barrier is present between wire and 
superconductor. Hence, here we study the regime of good 
electrical contact and equal potentials between S and the 
wire. 

We assume that transport is diffusive (scattering 
length being smaller than other length scales) so that 
the quasiclassical kinetic theory in the dirty limit can be 
applied. The Keldysh-Usadel diffusion equation for the 
Green's function, in Keldysh-Nambu space, in the ab- 
sence of time-dependencies, magnetic field and inelastic 
self-energy terms, can be written as 
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E is the energy measured from the chemical potential, 
and T) is the diffusion constant. The current is given by 
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(2) 



where Rn is the normal state resistance. 

The quasiclassical modeling of the diffusive transport 
through such a NSN structure can be splitted in solving 
the retarded and Keldysh parts of Eq. ([T]) respectively. 



A. The proximity effect 

In order to calculate the retarded part of the Green's 
function, G^, it is convenient to use the standard 9- 
parametrization, G^{x) = f3C0s6'(a;) H- f2 sin6'(a;). The 
function 9{x) is a measure of the superconducting corre- 
lations at a given point within the structure and satis- 
fies the Usadel equation^ 
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V—B{x) + 2iEsm0(x) = 0. 



(3) 



At the NS interfaces at a; = ±(i/2, the function 9{x) 
satisfies the following boundary conditional^ 
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= ±sm{0s-eN), x^±d/2, (4) 
= ^s^, ±d/2. (5) 



where S^n^s — \J Dj^^sjl-nT^ are the coherence lengths 
and Dat.s are the diffusion coefficients in N and S re- 
spectively. The proximity effect parameters 7 and 75 
are defined as 75 = Rb/pn^n and 7 = Ps^s/Pn^n, 
where Rb is the interface resistance and pn.s are the 
resistivities of the iV and S metals. These parameters 
have a simple physical meaning^^: 7 is a measure of the 
strength of the proximity effect between the S and N 
metals, whereas 75 describes the effect of the interface 
transparency. From here on, 7s ^ 1 is assumed, corre- 
sponding to a small barrier transparency. 

Solutions to the proximity effect problem in diffu- 
sive junctions have been extensively discussed in various 
regimes.^^'^^ Generally, a minigap exists in N, of the or- 
der of the Thouless energy. In the considered case of bulk 
N {dpj ^ ^jv), the minigap vanishes. The quasiparticle 
density of states is given by ReG = Re(cos^^) and the 
pair density of states is defined as ReF = Re(sin6'). It is 
straightforward to solve the Usadel equation ^ with the 
boundary conditions (U) and (O numerically. The calcu- 
lated densities of states RcGat and KcFn at the N side 
of an NS interface are shown in Fig. 2. It is seen that 
the quasiparticle spectrum in N is gapless while strong 
superconducting correlations exist at low E, described 
by ReF/v. As will be shown below, the existence of a 
nonzero Re F/v [E] at the NS interface is an essential in- 
gredient to our solution of the nonlocal conductance in a 
diffusive NSN structure. 



B. Distribution functions 

The Keldysh part of Eq. ([T]) provides the distribution 
of quasiparticles over energy. G^^ can be parametrized 
as G^ = G^f — fG^, where the distribution function 
/ can be split into parts that are, respectively, odd and 
even in energy, / = /l1 -I- /t^s- The kinetic equations 
for the longitudinal {/l) and transverse (fx) distribution 
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FIG. 2: Calculated density of states (ReGjv, dashed line) and 
pair amplitude (ReFjv, solid line) as induced by a supercon- 
ductor into a normal metal, with 7s=5, 7=0.1, and cIm 3> £,n- 
The functions are plotted for the JV-side of the interface, as 
indicated. 



functions have the forn 
V(Z?lV/l) 



Im/sV/i = 2/tA, (6) 
Im /sV/t - 0, (7) 



a.t X = ±d/2, where /Tiv(±<i/2) are the transverse distri- 
bution functions in the normal reservours. 

Note, that if the assumption 7 <C 1 is violated, the 
density of states Re G in the wire becomes finite at sub- 
gap energies. Then, the term Rei^ReFjv in the above 
boundary condition will be substituted by Re F Re Fn + 
Re G Re Gat. However, this will not lead to any qualita- 
tive changes in our results and the physical mechanism 
for the crossed Andreev transport remains the same in 
this case. 

Because of the small barrier transparency, the poten- 
tial mainly drops at the interfaces and we can assume 
the electric potential of the S'-reservoir and the wire to 
be zero, and the normal electrodes Ni^2 to be in equi- 
librium with the potentials Vi,2 respectively. Then, the 
distribution functions in A'^i 2 are 



fxiM = ^ tanh 
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2kBT 
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The kinetic equation has the solution fx = Ae^l^ + 
B e~^/^ w here the coherence length ^ is given by ^ = 
^JF)tJ2K^ which describes the conversion of quasiparti- 
cle current into supercurrent. This supercurrent is ex- 
tracted by the S'-reservoir at a; = 0. 



where Dt = (Re G)^ -I- (Re F)^ Dl = (RcG) - (ImF)^ 
and A is the gap inside the superconductor, /t and 
determine the quasiparticle and energy flow as can be 
derived^i from Eq. ([2]), giving respectively 
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The spectral supercurrent is given by Im/5 = 
iTr[f3(G^VG'" - G^VG^)] = Imi^^ReF^Vx, where 
X is the superconducting phase. We consider the regime 
when superconductivity in the S-wire is not influenced by 
the normal contacts A^i,2, which is realized when 7^1 
(large normal-state resistivity of N metal compared to 
that of S').23 Then, the product Imi^f Rei^J* in S is 
nonzero only at a narrow energy range near A. Since we 
are interested only in the energy range £' < A , this al- 
lows us to neglect the terms with Im Ig cx Im F^ Re Fg- 
in the kinetic equations ^ and ([7]) at these energies, and 
the equations for Jt and /i decouple. The equation for 
Jt in the wire, consequently, becomes 



= 2/tA, 



(10) 



with the boundary conditions at the NS interfaces given 
bj^ 

DTjB^fT = ±'ReFReFNifT-fTm) (11) 



III. RESULTS AND DISCUSSION 

Solving Eq. (jlOp with the boundary conditions, Eq. 
(fTTj) . the solution of fr in the case of two symmetric NS 
interfaces is given by 
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B = N 



hie-"/^^ (7B -N) + fT2e''/^^ (7b + N) 
e-d/i (7B - Nf - e'i/i (75 + Nf ' 

/Tie'^/^^ (7B + N) + fTie""/^^ (75 " N) 
e-d/i _ Nf _ grf/? (75 + Nf 



(12) 



where TV = Re Re F^r . 

The bias condition provides the final equation from 
which the distribution functions and currents are derived. 
The first interface is biased with a voltage Vbias, from 
which /ti is then known. Two bias conditions for the 
second interface are considered, (a) Iqp — 0: For a zero 
total current through the second interface, i.e. an open 
connection, or an ideal Voltmeter, the nonlocally induced 
voltage has to be found self-consistently from Eq. (E 
and the additional boundary condition. 



DT-^frdE 
ox 



0. 



(b) Vs = Vn2- Under the alternative bias condition 
of zero potential difference across the second interface, 
/t2 — and Eq. p2|) directly provides the solution for 
/t in the superconductor. 

Note, that it is essential to our approach that ReFjy 
is nonzero at subgap energies. When these proximity 
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FIG. 3: Spectral quasiparticle current at zero temperature 
accross the second interface under a voltage bias of the first 
interface, Vbias = 0.3A, and a zero total current across the 
second interface as sketched. The superconducting interlayer 
thickness d = 0.3^, and 7s = 5. Inset: The response of the 
induced nonlocal voltage across the second interface (which 
is the measured quantity of Ref,9) as function of the applied 
voltage bias across the first interface. 



FIG. 4: Nonlocal spectral quasiparticle current at zero tem- 
perature across the second interface under a voltage bias of 
the first interface, Vbias = 0.3A, and a zero voltage bias of 
the second interface as sketched. 7s = 5 and d — 0.3^ (solid 
line), d — ^ (dashed line), d = 3^ (dotted line). Inset: The 
response of the total nonlocal current across the second in- 
terface as function of the applied voltage bias across the first 
interface, for 73 = 5 and d — 0.3^. 



induced correlations are neglected, the quasiparticle cur- 
rent Iqp and the nonlocal effects vanish, coinciding with 
the results from previous tunnel Hamiltonia n^^i^^ and cir- 
cuit theorjsi^ models. 

The resulting spectral quasiparticle currents across the 
second interface are shown in Figs. [3] and 01 for a thin 
superconductor of c? = 0.3^ embedded symmetrically be- 
tween two tunnel barriers with 7_b = 5. In the case of 
bias condition (a), a nonlocal voltage, Vni, is induced in 
the unbiased normal metal electrode, see Fig. [31 as was 
experimentally observed.^ Vni is the source of a local An- 
dreev reflection process at the second interface, as char- 
acterized by the low-energy peak in Iqp (E). In Fig. [U 
also Iqp (E) across the first interface is shown, to illus- 
trate the different quasiparticle tunneling and reflection 
processes. A negative current corresponds to a flow of 
electrons in the positive direction. Iqp (E) at the second 
interface, outside the spectral region of Andreev reflec- 
tion ([£'1 > Vni) is physically caused by the two non- 
local processes of elastic cotunneling (EC) and crossed 
Andreev reflection (CAR). The reverse backflow process 
of EC and CAR from the right to the left electrode re- 
sults in a suppression of Iqp (E) at the first interface for 

\E\ < Vni. 

In bias situation (b), only the nonlocal currents are 
contained in Iqp {E), see Fig. |4l The magnitude of the 
nonlocal currents scales exponentially with d/^ [see so- 
lution Eq. and Fig. 2], as expected from the fact 
that a Cooper pair has size ^ and that nonlocal effects 
exist by virtue of the coupling to the superconducting 
condensate. Thus, the Thouless scale is not relevant to 
our model. The nonlocal currents and voltage scale with 
7^^, which provides a much larger effect than was found 



from the tunnel Hamiltonian approach (fourth power in 
transparency) 

In the linear response regime (low temperature and low 
voltage), it was shown recently by Morten et alr^ that 

— ^ = Geo (el4ias) ± Gcar (el4ias) , (13) 

O Vbias 

where Geo and Gcar are defined as the nonlocal con- 
ductances for EC and CAR respectively. From the fact, 
that Z?L = in the superconductor at subgap volt- 
age, and Eq. it follows that II — 0, and Gec — 
—Gcar- This means that the nonlocal quasiparticle cur- 
rent, Iqp {E) = - (Gec - Gcar) /ti, as shown in Fig. gl 
is carried by the EC and CAR processes equally, i.e. 
Iec = ^AR = Iqp (E) /2. The CAR hole in the sec- 
ond electrode can, therefore, be thought of as moving in 
the negative direction, while the EC electron moves in 
the positive direction (see Fig. 1). 

A treatment of the nonlocal processes in terms of elec- 
trons and holes can be derived from the respective elec- 
tron and hole distribution functions, /e and fh, that 
are given by2i ^ = (1 ± - fr) /2, when also /l is 
treated self-consistently. 

The sign and magnitudes of the modeled nonlo- 
cal effects are of use for the interpretation of recent 
experiments;2iiiii2 although many aspects of the exper- 
iment related to the geometry are not covered by our 
model. The obtained equal contribution of EC and CAR 
to nonlocal currents in NSN structures indicates that 
additional quasiparticle manipulations are necessary be- 
fore the device can be considered as a useful source of 
entangled particles. Creating a non-equilibrium distri- 
bution in the electrodes (for example by energy-filtering 
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in a Fabry-Perot structure) in this respect would be beni- 
ficial. 



parency, providing large nonlocal effects, of the order of 
experimentally observed magnitude. 



IV. CONCLUSION 

In summary, within the assumptions of quasiclassical 
Keldysh-Usadel theory, we find that nonlocal EC and 
CAR effects exist in a diffusive quasi- ID wire by virtue 
of the proximity effect. We have found that CAR and 
EC have the same sign in their contribution to the non- 
local spectral quasiparticle current and nonlocal voltage. 
CAR and EC scale with the square of the barrier trans- 
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